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Abstract 

We consider a Friedmann brane moving in a bulk impregnated by radiation. The setup is strongly 
asymmetric, with only one black hole in the bulk. The radiation emitted by this bulk black hole can 
be reflected, absorbed or transmitted through the brane. Radiation pressure accelerates the brane, 
behaving as dark energy. Absorption however generates a competing effect: the brane becomes 
heavier and gravitational attraction increases. We analyse the model numerically, assuming a total 
absorbtion on the brane for k = 1. We conclude that due to the two competing effects, in this 
asymmetric scenario the Hawking radiation from the bulk black hole is not able to change the 
recollapsing fate of this brane-world universe. We show that for light branes and early times the 
radiation pressure is the dominant effect. In contrast, for heavy branes the self-gravity of the 
absorbed radiation is a much stronger effect. We find the critical value of the initial energy density 
for which these two effects roughly cancel each other. 
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I. INTRODUCTION 



According to brane-world models, our universe is a hypersurface (the brane) embedded 
in a 5-dimensional spacetime (the bulk), in which Einstein's gravity holds. The attempts 
to increase the number of dimensions in physics have a long history, originating with the 
pioneering works of Kaluza and Klein. In the more recent models presented in [1], the (flat) 
extra dimensions are still compactified. The novelty of the brane-world models generalizing 
the original proposal [2] of Randall and Sundrum (RS) consist in allowing for one additional 
noncompact dimension. Instead of compactifying, the bulk is warped. 

Gravitational dynamics on the brane is given by a modified Einstein equation, which 
arises from the projection of the 5-dimensional dynamics and the junction conditions across 
the brane. A major difficulty in this theory is that due to a new source term, arising from 
the projection of the bulk Weyl tensor to the brane, the system of differential equations 
describing gravitational evolution is not closed on the brane. This new gravitational theory 
allows general relativity as the low energy density limit (as compared to the brane tension 
A). 

Brane-world models were initially Z 2 -symmetric, with identical copies of the bulk on the 
two sides of the brane. Equivalently, the brane could be imagined as a domain boundary. 
This assumption, although it considerably simplifies the model, seems unnecessarily restric- 
tive in the context of brane-world models with curved branes (generalized RS type 2 models). 
In this context, our observable universe may be imagined as a Friedmann brane moving ei- 
ther in a Schwarzschild-anti de Sitter bulk (if the bulk is vacuum), or in a 5-dimensional 
Vaidya-anti de Sitter (VAdS5) spacetime (if the bulk contains radiation). The modified 
Einstein equation in the ^-symmetric case was given covariantly in [3], and generalized 
for asymmetric embedding in [4]. The effects of asymmetric embedding were also analyzed 
in a more generic class of models, containing induced gravity contributions [5]. There it 
was shown that the late time universe behaves differently under the introduction of a small 
asymmetry of the embedding: on the two branches of such theories. While on the RS branch 
asymmetry produces late-time acceleration, on the DGP branch the late-time acceleration 
is diminished by asymmetry. 

Whenever there is radiation in the bulk, it can arise from both brane and bulk sources. 
The situation when the brane radiates into the bulk was already considered both in a 
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symmetric [6] and in an asymmetric scenario [7], [8] in the framework of generalized RS 
theories. In both scenarios a radial emission of radiation into the bulk was considered. 
A more realistic set-up, allowing for the emitted gravitons to follow geodesic paths, was 
advanced recently in [9]. 

In this paper we discuss the other possibility, when the radiation is emitted by sources in 
the bulk. For this, we embed the Friedmann brane asymmetrically into VAdS5 spacetime. 
The symmetry of the embedding is severely broken by allowing for a radiating black hole on 
one side and no black hole on the other side of the brane. Due to this asymmetry, there is 
a radiation pressure acting on one side of the brane. This radiation pressure continuously 
accelerates the motion of the brane, leading to accelerated cosmic expansion. In principle, 
this mechanism can be a dark energy candidate. 

In brane-world models, standard matter is confined to the brane. Therefore the radiation 
coming from the bulk should be of non-standard nature. We discuss here the case when 
this is the Hawking radiation of the single bulk black hole. The expression of the Hawking 
radiation in a 5-dimensional Schwarzschild-anti de Sitter (SAdS5) spacetime was derived for 
the curvature index k — 1 in [10]- [12]. This energy density was employed, but for k = in 
[13] for the study of asymmetrically embedded Friedmann branes into VAdS5. 

Some part of the radiation reaching the brane will be absorbed, some part will be re- 
flected and the remaining part will go through. We will study here only the models without 
reflection in order to have the VAdS5 spacetime in both bulk regions (Fig 1). We disregard 
the possibility of reflection, because there is no known exact solution with cosmological con- 
stant describing a crossflow of radiation streams. Such a spacetime with both incoming and 
outgoing null dust streams, but without a cosmological constant is known in 4-dimensions 
[14], however no generalization including a cosmological constant has been found. Thus we 
discuss the case when the reflected component can be neglected. 

We keep however the absorption, as an essential element of our model, in contrast with 
the one presented in [13], where there is a black hole on each side of the brane, but their 
radiation is completely transmitted through the brane. 

By retaining the absorption on the brane, we obtain novel features. The absorbed ra- 
diation appears on the brane as continuously emerging energy. Thus the energy density 
on the brane increases, strengthening the gravitational self-attraction of our universe. As 
consequence of absorption, cosmic expansion slows down. 
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Therefore there are two competing effects in our model: radiation pressure accelerates, 
while the absorbed radiation decelerates cosmic expansion. We would like to study how 
their balance affects the Friedmann brane-world. We derive the relevant equations of the 
model in Section II and we comment on the physical interpretation in Section III. 

As our model emphasizes the effects of the absorption (as opposed to [13]), we finally 
choose a total absorption on the brane in Section IV. We also employ k — 1, as the formula 
for the energy density of the Hawking radiation was derived for this case [10]- [12]. Then we 
introduce new dimensionless variables, adapted to our choice of k — 1. 

We provide quantitative results by numerical analysis in Section V. First we discuss 
cosmological evolution in the case of a non-radiating bulk black hole. Then, by switching on 
the radiation, we see how the two competing effects modulate the cosmological evolution. 
We show that a critical behaviour can appear, when these two competing effects roughly 
cancel each other. We discuss the relevance of these results in a broader context in Section 
VI. 

Throughout the paper a tilde distinguishes the quantities defined on the 5-dimensional 
spacetime and a hat denotes dimensionless quantities and units with c—l — h are chosen. 



II. FRIEDMANN BRANE IN VAIDYA-ANTI DE SITTER BULK 



The Friedmann brane representing our observable universe is embedded in the VAdS5 
spacetime: 



where 



and 



ds 2 = —f (v, r; k) alv 2 — 2dvdr 

+ r 2 [d X 2 + H 2 (x;k){d6 2 + sin 2 Odcf 2 )] , 



(1) 



n( X ; k) 



sin x , k = l 
X, k = 
sinh x , k = — 1 



(2) 



/ ( v ,r; k)=k- 



2m (v) k 2 A 
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(3) 



Here k 2 = 87rG( 5 ) is the gravitational coupling constant in the bulk (G(5) representing the 
5-dimensional gravitational constant) and it has the dimension of (length 4 mass _1 time~ 2 ). 



The null coordinate v is ingoing (the v = const, lines are outgoing). As a simplifying 
assumption and because we would like to focus on the radiation effects, we choose the same 
bulk cosmological constant A on the two sides of the brane. The mass function m (v) is freely- 
specifiable, however different on the left and right sides of the brane: (v) ^ (v). We 
note that, as defined here, the mass function m(v) has the dimension of length 2 . This is, 
because it denotes G^M (v) /c 2 , with M (v) the true mass. Since we have chosen h = 1 = c, 
dimensionally speaking, length, time and mass -1 are the same. Thus technically the mass 
function m (v) represents the cumulated mass of the bulk black hole and radiation multiplied 
by G(5) (of dimension mass" 3 ). 

The source of the metric (1) is pure radiation considered in the geometrical optics limit 
(null dust), with energy- momentum tensor 

f» D = i/;(v,r)l a l b . (4) 
Here ip (v, r) is the energy density and / is a (future-oriented) null 1-form: 

/ = — — — n + u . (5) 

v 

Here n is the dual 1-form of the unit normal to the brane and u the dual 1-form of the 
timelike (negative) unit vector. We define these vectors in what follows. We give the brane 
by the embedding relations v — v (r) and r = a (r), with r the cosmological time. Then we 
choose the timelike unit vector u such that r is adapted to u: 

d . d . d . . 

<9r dv dr 

(A dot denotes derivatives with respect to r). As consequence of u a u a = —1, we have on 
the brane 

fv = -a+(a 2 + f) 1/2 . (7) 

Here we have chosen the positive root in order to have v positive. Now we can write the 
1-form field u a as 

u = - (r 2 + f) 1 ' 2 dv - vdr = -dr . (8) 
The unit normal 1-form to the brane can be expressed as 

n = —rdv + vdr . (9) 
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FIG. 1: The radiation from the left bulk black hole is partially absorbed on the brane and partially 
transmitted through into the other bulk region containing no black hole. 

The bulk Einstein equation for the metric (1) and the source term 

fL cd = -Ag ab + f% D , (10) 

establishes a connection between the energy density if) and the time derivative of the mass 
function m = (dm/dv) v as 

where v is given by Eq. (7). 

The square root from v can be eliminated by using the relations derived in [4]: 

(a 2 + fn, L f 2 = ±B AB 
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B = — ^0 + A ) > 



6Am , . 

= v^TT) ' < 12) 

(In the above and the forthcoming formulae the first subscript refers to the upper sign and 
A always is understood as the difference taken between the R and L regions. An overbar 
denotes the average of a physical quantity, taken on the two sides of the brane.) Therefore 
Eqs. (7), (11) and (12) imply 



Wr,l — — ~ 2 on — 
« • a Jr,l 



k a , , . 3 Am 

-a T —r- (p + A) + 



6 vr ' n 2 a 3 (p + X) 



(13) 



Why do we need such a detailed expression for ipR^t The answer is that radiation energy 
density appears in two of the key equations governing cosmological evolution. (These equa- 
tions emerge from the modified Einstein equation.) The energy densities ipR t L are present 
both in the energy-balance equation 

p + 3-(p + p) = -Aij , (14) 
a 



and in the Raychaudhuri equation: 

a Aq k 2 
a 3 6 



P 



1 + t) + *( 1 + 9 



2m 27 (p - A) (Am) 2 
'~^ + K 4 a 8 (p + \f 
k 2 ^ 3AmAip 



3 Pa 4 (p + A) 2 ' 

However, the third essential equation, the Friedmann equation is independent of tp: 

a 2 + k Aq | n 2 p / | p 



(15) 



a 2 3 ' 3 V" ' 2A 

2m 9 (Am) 2 . . 

+ (16) 

(The above equations are particular cases of the equations derived in [4] in a more generic 
set-up.) In the above two equations A denotes the 4-dimensional cosmological constant 
given in terms of the brane tension A and bulk cosmological constant as: 

2A = k 2 A + £ 2 A . (17) 

The system of equations (13)-(16) governs the motion of the brane, seen from the brane 
point of view as cosmological evolution. 

III. THE FREE FUNCTIONS OF THE MODEL 

There are a number of free functions in our model. Among them, by straightforward 
algebra on Eq. (13) the quantities Aip and ip can be expressed in terms of m, Am, m = m 
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and (Am) = Am as 



2k 2 a 3 
1 



3Am 



a — 



K 2 a 3 (p + A) 
(p + A) (2mF + + AmF_) , 



(2mF_ + AmF+) 



^ = 



4/t 2 a 3 
1 



3Am 



ft- 



! a 3 (p + A)_ 



(2mF + + AmF_) 



8a 2 



(p + A) (2mF_ + AmF + ) , 



where we have introduced the shorthand notations 

1 1 _ 2a 2 Am 
f R h ~ [a 2 ]) 2 - (Am) 
1 1 2a 4 J 



a 2 / 



fR h (a 2 f) 

K 2 \ A 

~6 3~ 



{Amf 



a + ka — 2m . 



(18) 
(19) 

(20) 
(21) 
(22) 



Then in Eqs. (14)-(16) the role of the free functions ip and Atp is taken by m and Am . 

In what follows, we discuss the physical interpretation of these free functions. First we 
show that Am is determined by three quantities. These are the energy density of the brane, 
the scale factor and Hubble parameter. Second, that ipL is given by the energy emission of 
the bulk black-hole, evaluated at the left side of the brane. Finally we interpret Am as the 
energy absorption on the brane. 



A. The energy content of the brane 



By rearranging the Friedmann equation, we find Am as function of the cosmological 
perfect fluid energy density p, scale factor a and Hubble parameter H = a/a: 



(Amf 



2k 2 
~3\ 



a 8 (p + A) 



(23) 



In the derivation we have employed the relation among the brane tension and the gravita- 
tional coupling constants in 5 and 4-dimensions: 



K 

T 



K 
If 



(24) 



It is easy to argue, that m# > m^ as m# contains in addition to m^ the contribution from 
the brane energy momentum. Therefore Am > 0. 



We stress here that Am is not strictly the mass of the brane. Rather, it represents the 
difference of the mass functions of the VAdS5 metric, evaluated on the two sides of the 
brane. Although, dimensionally speaking, Am is mass -2 , it is directly proportional to the 
difference in the cumulated masses M (v) of the bulk black hole and radiation on the two 
sides of the brane. The reason why jjir ^ uil is that the brane contributes toward the mass 
function both though its energy momentum and its curvature. In a less strict sense then 
Am can be regarded as a measure of the energy content of the brane. 



B. Hawking radiation in the bulk 

The rate of decrease in mass of an evaporating black hole in an SAdS5 spacetime with 
k = 1 was computed in [10]- [12]. In [13] the assumption was advanced that this agrees with 
the rate of decrease of the mass function m^ (v) in the VAdS5 spacetime, even for arbitrary 
k. Numerical results were given in [13] for k — 0. We adopt here the same identification, 
but only for k — 1, for which the proof of [10]-[12] holds. Then the energy density of the 
radiation escaping from the black hole, evaluated at the left side of the brane is: 

, L = £(5) , (25) 



24ir 9 a 3 mL 



d + (a 2 + / L ) 1/2 

where ( is the Riemann-zeta function. Here m^ contains both the mass of the black hole and 
the energy of the Hawking radiation. This radiation spreads away with the velocity of light 
and it overtakes the brane, which has a subluminal motion, therefore rh^L < 0. Moreover, 
as part of the energy radiated away is captured on the brane, there is less radiation escaping 
in the R region than there would be in the absence of the brane. Therefore tHl < m R < 
and in consequence Am > 0. For the same reason Aip < holds. 



C. Energy absorption on the brane 

The decomposition of the pure radiation energy-momentum tensor (4) with respect to 
the brane can be done by employing Eq. (5). We obtain: 

(Tr,E) = i>R,L {u a u b + 2u {a n h) + n a n b ) . (26) 

V / ab 

On the L region the radiation escaping the bulk black hole is characterized by ip L . Some of 
this radiation will pass unaffected through the brane into the R region (V\r)- Then we have 

9 



to interpret the difference of the transmitted and incident energy-momentum tensors: 

AT a ^ D = Atp (u a u b + 2u (a n b) + n a n b ) 
= - {p rad u a u b + 2q\fu b) + p rad n a n b ) 5 (y) . (27) 

Therefore the difference in the energy density of the radiation arriving from left and escaping 
into the right regions appears on the brane as a dust with energy density p rad ~ (— A^) > 0. 
The energy flux ariving from the bulk into the brane is q r a ad ~ (—Aip) n a . Finally, the 
radiation pressure p rad ~ (—Atp) > drives the brane into an accelerated motion. This in 
turn accelerates the cosmic expansion. 

With A^ < 0, the acceleration caused by the radiation pressure is manifest from the last 
term of the Raychaudhuri equation (15). 



IV. TOTAL ABSORPTION ON THE BRANE 



In this section we consider a simple model, in which k = 1 and A = 0. We allow the 
brane to absorb all of the radiation escaping from the bulk black hole. With this choice the 
equations simplify considerably and we can study the effects of the absorption, which were 
not taken into account in previous treatments. 



A. Brane dynamics 



For total absorption ipn = and rriR becomes a constant. Two free functions are then 
left, which can be chosen as m := Am and ip := ipi- The energy-balance and Raychaudhuri 
equations (14) and (15) become 





a 
a 



(28) 



K 



p 



2p 
A 



2rriR m 27 (p — A) (m) 



a 4 



+ 



a 4 W (p + XY 
3ipm 



(29) 



6 ' K 2 a 4 (p + \y 

We remark that with the assumption ip R = both ip and Aip become determined. The 
Friedmann equation (16) emerges now as a consequence of the energy-balance equation 
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(28) and the Raychaudhuri equation (29). Stated otherwise, only two of the Friedmann, 
Raychaudhuri and energy-balance equations are now independent. (The interdependence of 
these three equations is a generic feature of the standard cosmological models, but not of 
brane-world universes. Indeed, prior to assuming the total absorption, the three equations 
were independent.) 

From Eqs. (7), (11) and (25) the variation in time of the brane mass function is found: 

5C(5)£ 2 . , 

This equation is an ordinary differential equation which can be integrated for t > t . The 
result is: 

m = m R - \J (m R - m ) 2 - 5 ^| 9 — (* ~ *o) , (31) 

where by m we have denoted the brane mass function before the Hawking radiation reaches 
the brane (this happens at t — to). The brane mass function m is monotonically increasing in 
time until tc = to + 36n 9 {tur — m ) 2 /5£ (5) k, 2 , when it reaches its maximal value m = uir. 
When this happens, the bulk black hole has already completely evaporated and the resulting 
Hawking radiation is entirely absorbed by the brane. 

The square root in Eq. (25) can be eliminated by employing Eq. (12). We obtain 

* = ■ (32) 



2W (m R - m) [i + (p + A) + ^J^] 



With this both free functions m and ip are given explicitly as functions of time. This has 
been achieved by specifying both the rate of absorption and the energy density of the black 
hole radiation Eq. (25). 



B. Dimensionless variables 

We introduce dimensionless variables in a slightly different way than in [6] and [7]. The 
reason for this is that the chosen comoving radial coordinate x nas to be dimensionless in 
the k — ±1 cases. Therefore the scale factor has to carry the dimension, unless in the case 
k — 0, discussed in [6] and [7]. The dimensionless variables introduced below are well suited 
for k — ±1 but they also apply for k = if in this latter case the scale factor, rather then 
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the x-coordinate is regarded dimensional. These variables are: 

jj 

t = Ct,a = Ca,H = — , 

P - V r i> 
p= A' P= A'^=CA' 

m = C 2 m , m R = C 2 m R . (33) 

Here C = k\/\ and its inverse represents a distance scale. (The dimensions of k 2 and A are 
length/mass and mass/length 3 , respectively.) We note here that the dimension of tp given 
by Eqs. (11) or (25) is consistent only by choosing h = 1 = c, implying that the dimension 
of mass is length -1 . Then C can be thought as a mass scale as well. Alternatively, by the 
choice c = 1 (therefore length and time having the same dimension), C becomes a basic 
energy scale. 

Written in the variables (33), the system of equations describing cosmological evolution 
contains the following dimensionless parameter: k s X 3 = 36/t 4 A = 23047r 2 A/M 4 , where M p 
denotes the 4-dimensional Planck mass. By denoting the brane tension as A = M 4 ,, the 
free parameter can be chosen as Mt/M p . However for later notational convenience we 
introduce the dimensionless parameter (3 = 40C (5) (M T /M P ) A / (3n 7 ). The energy-balance 
and Raychaudhuri equations (14) and (15), written in terms of the dimensionless variables 
(33) and parameter f3 read 

= p + 3H(p + p) -f , (34) 
H' = -H 2 - \p{\ + 2p) - \p{\ + p) 



2m R in 9(p — l)(m) 



2 



a : o 4 2a 8 (p + if 
V6a 4 (p + lf 



^ + fT^A — 7^2 > ( 35 ) 



with 



and 



# = rJ—, 77— r (36) 



a 4 {rfiR — m) 



tt , p+l , 3m 



V6 {m R - m) 

Here a prime denotes derivates with respect to the dimensionless time t . Eq. (37) integrated 
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gives the temporal evolution of the dimensionless mass function of the brane 



m = m R -^J (m R - m ) 2 - -^L (?- f ) • (38) 

The generalized Friedmann equation (16) in dimensionless variables readily follows: 

^9 1 P / p\ 2m R 



a 2 3 V 2 J a 4 

m 3 (m) 2 
-~ + — — n • 39 

a 4 2a«(p+l) 2 1 ; 

V. NUMERICAL RESULTS 

In this section we give numerical results on the evolution of the brane-world universe in 
the presence of one radiating bulk black hole. We assume the movement of the brane starts 
on the apparent horizon 

r AH = + ^~+\2m L (40) 

of the bulk black hole (therefore t — 0). This horizon is dynamical, as m L decreases due 
to the escaping Hawking radiation. The brane leaves the horizon at m L0 = m R — m . In 
order to have clear graphical solutions, we have found convenient to chose the dimensionless 
mass values of m R = 0.3 and mo = 0.29. Then the initial position of the apparent horizon 
is at 0.1411870178. Away from the horizon, the brane absorbes the Hawking radiation. We 
investigate a radiation dominated universe, so p = jo/ 3. We choose the free parameter as 
M T /M P = 1/10. In what follows, we compare the evolutions of the brane in the absence 
and in the presence of the radiation coming from the bulk black hole. 



A. The evolution of the brane-world without bulk black hole radiation 

In the simplest case when we switch off the radiation from the bulk black hole, beyond 
the standard density and curvature source terms in the Friedmann equation (39), the model 
includes dark radiation and a strong asymmetry due to the chosen setup of a single left 
black hole in the bulk. We would like to see the evolution of such a brane universe, for later 
comparison with the radiating case. With the Hawking radiation from the bulk black hole 
switched off, in = const. 
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FIG. 2: The evolution of the scalefactor in the absence of the Hawking radiation shows the char- 
acteristic pattern of a closed universe. The plot is for the initial brane density po = 520. 




FIG. 3: Evolution in the absence of the Hawking radiation of the Hubble parameter (dotted line) 
and brane energy density (solid line). The initial brane density is po = 520. 

Our numerical analysis shows that none of the brane effects (quadratic energy density, 
dark radiation and asymmetry source terms) will change significantly the usual behaviour of 
a closed (k = 1) universe, as can be seen from the behaviour of the scale-factor, represented 
in Fig 2. The Hubble parameter and the brane energy density evolve accordingly. The 
Hubble parameter continuously decreases, while the energy density tends to infinity both at 
the beginning and at the end of this brane-world universe. These behaviours are shown in 
Fig 3. 
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FIG. 4: Evolution of the dark radiation (DRST, dotted line) and asymmetry (AST, solid line) 
source terms of the Friedmann equation (for the initial brane density po = 520). 



The evolutions of the dark radiation source term 

2uir rh 



a 4 a 4 ' 



(41) 



and of the asymmetry source term 



3(a)2 2 (42) 
2a 8 (p + l) 2 

are represented in Fig 4. Both contributions are positive, so they act in a similar way to 
ordinary matter sources in the Friedmann equation. The dark radiation dominates over the 
asymmetry during the whole evolution, the difference being less when the universe reaches 
its maximal size. As the initial value of the brane energy density is lowered, these two source 
terms become comparable (Fig 5). For even lower initial energy densities the asymmetry 
source term dominates at both early and late times (Fig 6). 



B. The energy density of the Hawking radiation 

The energy density if) (evaluated near the brane, on the left region) of the Hawking 
radiation from the bulk black hole evolves; cf Fig 7. In the expanding era ip shows a sharp 
decrease. This is due to the fact that the distance between the brane and the horizon 
increases. Interestingly ip starts to increase before the maximal size of the closed universe 
is reached. Later on ip keeps increasing as the universe contracts. The higher values of 
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FIG. 5: Similar evolutions of the dark radiation (dotted line) and asymmetry (solid line) source 
terms in the Friedmann equation at po = 250. 
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FIG. 6: For light branes (j>q = 100) the asymmetry (solid line) source term dominates over the dark 
radiation (dotted line) source term in the Friedmann equation for both early and late universes. 

ip towards the end of contraction mathematically can be explained by the presence of the 
forever decreasing factor H in the denominator of ip (see Eq. (36)). Stated otherwise, from 
Eq. (11) ip~v so that from a physical point of view, the increase in ip at late times, exceeding 
its value at early times can be explained as a Doppler effect. The contracting brane absorbs 
more energy in unit time than the expanding brane at the same scale factor value. 

We have checked, that with the Hawking radiation switched on, the basic behaviour of 
the closed universe will be very similar to that in the case without this radiation, discussed 
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FIG. 7: The evolution of the energy density of Hawking radiation for po = 520. 

in the previous subsection. Even in the presence of the radiation the cosmological evolution 
shows the same patters. Although there are modifications, which will be discussed in detail 
in the following subsection, they behave as small perturbations during the whole evolution 
of this asymmetric brane-world universe. 

C. The evolution of the brane-world with the bulk black hole radiation 

The evolution of the Hubble parameter, the scalefactor, the brane energy density and 
different source terms of the Friedmann equation, in the presence of the radiation from the 
bulk black hole is very similar to their evolutions in the absence of the radiation. This is due 
to a delicate balance between two effects. First, the brane is pushed away by the radiation 
pressure, which tends to accelerate it. Second, as all of the radiation is absorbed on the 
brane, the self-gravity of the brane increases, which tends to decelerate it. Still, there are 
interesting, through small effects due to the radiation. We illustrate them by representing 
the differences in the scalefactor, the brane energy density and different source terms of the 
Friedmann equation in the radiating and non-radiating cases. More specifically, we denote 
by Af the quantity A in the radiating case and we plot the evolution of S A = A^ — A . 

In Fig 8 we show the evolution of the differences in the scalefactors when the radiation is 
switched on and off, for three different initial brane energy densities. Surprisingly, we have 
found a critical-like behaviour for p = 520. The difference in the scalefactors is extremely 
close to zero, but it shows a sinusoidal-like pattern. For higher initial densities the scalefactor 
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FIG. 8: The difference between the scale factors in the radiating and non-radiating case versus time. 
For po = 520 a critical-like behavior is observed (solid line) , when the increase of the self-gravity of 
the brane due to absorbtion is roughtly compensated by the radiation pressure. For initial densities 
much smaller (e.g. po = 100, upper dotted line), the effect of the radiation pressure is dominant. 
For higher initial densities than the critical one (po = 2000, lower dotted line) the increase of 
self-gravity due to absorbtion dominates over the radiation pressure. 

is smaller in the radiating case throughout the cosmological evolution. This means that the 
self-gravity of the absorbed radiation dominates over the radiation pressure. For initial 
pressures higher than the critical one, the scalefactor is higher in the radiating case during 
the whole cosmological evolution. This means that the radiation pressure dominates over 
the self-gravity of the absorbed radiation. 

The behaviour of the difference in the energy densities in illustrated in Fig 9. Among 
the previous graphs, Figs 2, 3, 5 and 9 were plotted for the critical initial density. For other 
density values these graphs would be merely scaled, and would show similar features. 

This is not true for the source terms of the Friedmann equation. We have seen various 
behaviours even in the non-radiating case (Figs 4-6), when the initial brane energy density 
is varied. In Figs 10-12 we have plotted the evolution of the difference, in the radiating 
and non-radiating cases, of the asymmetry and dark radiation source terms, for the critical 
initial energy density (Fig 10), for a much lighter brane (Fig 11) and for a much heavier 
brane (Fig 12). With the exceptions of the very early and very late stages of the evolution 
of the universe, the two contributions roughly cancel each other for the critical initial energy 
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FIG. 9: The difference between the brane energy densities in the radiating and non-radiating case 
for po = 520. 
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FIG. 10: The difference between source terms in the Friedmann equation computed when the 
radiation is switches on and off. The asymmetry source term (5 AST, solid line) is increased, while 
the dark radiation source term (5DRST, dotted line) is decreased by radiation. The plot is for the 
"critical" initial density po = 520. 

density. For light branes the increase in the asymmetry term at early times is greater than 
the decrease in the dark radiation term. For heavy branes, by contrast, the decrease in the 
dark radiation term due to the Hawking radiation from the bulk exceeds the increase of the 
asymmetry term. 
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FIG. 11: The same as in Fig 10, but for po = 100. The main difference with respect to Fig 10 
is that at early times the increase in the asymmetry term is faster than the decrease in the dark 
radiation term. 
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FIG. 12: The same as in Fig 10, but for po = 2000. The decrease in the dark radiation term 
dominates over the increase of the asymmetry term. This effect is more accentuated than at the 
critical initial density, shown on Fig 10. 

VI. DISCUSSION 

Table-top experiments [15] on possible deviations from Newton's law currently probe 
gravity at sub-millimeter scales and as a result they constrain the characteristic curvature 
scale of the bulk to I < 0.1 mm. Expressed in units of inverse energy (when c = 1 = K) this 
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curvature scale is 

/max = 506.7 7 eV" 1 = 0.5067 7 x 10 12 GeV" 1 . (43) 

This is usually known as / < 10 12 GeV -1 [16]. 

The 4-dimensional coupling constant k 2 and 4-dimensional gravitational constant G are 
related to the 4-dimensional Planck mass M P as k 2 = 8nG = 8n/Mp, with M P ps 10 19 GeV. 
The 5-dimensional Planck mass is defined as k 2 = 8n/M? 5 y As M( 5 ) depends on both M P 
and on the characteristic curvature scale I as Mh, = Mp/l [2], we get 

M (5) min = 6. 65 x 10 8 GeV . (44) 

This result is referred usually as M P > M( 5 ) > 10 8 GeV [16]. Due to the presence of the extra 
dimension accessible via gravity, brane-world theories allow for the effective 4-dimensional 
Planck scale on the brane to be be much higher than the true, 5-dimensional Planck scale. 

The knowledge of both / and M P gives a characteristic lower limit for the brane tension 
in the two- brane model [17]. In units c—l — h this is 

"\ M 2 

Ami„ = = 138.59 TeV 4 . (45) 

This result is frequently quoted as A > 1 TeV 4 [16]. On the other hand Big Bang Nucle- 
osynthesis (BBN) constraints result in a much milder lower limit, A > 1 MeV 4 [18]. There 
is also an astrophysical limit for A, this being sensitive to the equation of state of a neutron 
star on the brane [19]. For a typical neutron star (with uniform density) this is A > 5 x 10 8 
MeV 4 , in between the two previous lower limits. 

Our choice for Mt/M p = 1/10 employed in the numerical analysis of the previous Section 
and the numerical value of the Planck-mass M P w 10 19 GeV implies M T ps 10 15 TeV. In 
consequence A = (Mr) 4 = 10 60 TeV 4 , well above all previously enlisted lower bounds. This 
choice is in agreement with the range of A represented in Fig. (1) of Ref. [18]. The reason 
we have chosen such a high value of the brane tension is the following. Decreasing A would 
mean to decrease Mt accordingly. As ip is proportional to (Mt/M p ) 4 (cf Eq. (36)), by 
choosing a significantly smaller value of A the effect of the radiation energy density becomes 
very small, in the range of the numerical errors. 

In this paper our intention was to focus on the effects of the asymmetric setup and bulk 
radiation absorbed on the brane. Therefore we have simplified the model by setting A = 0. 
By means of our choices of m R and m the dark radiation term was chosen of comparable 
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size to the matter energy density source term, slightly above the limits established in [20] 
and [21] from BBN for the dark radiation term. We note however that those limits were 
derived in a symmetric setup, and in the presence of a cosmological constant, which is not 
the case here. 

The question arises of whether the dark radiation term can take the role of a cosmological 
constant producing cosmic acceleration. Contrarily to optimistic expectations, this is not 
exactly happening in brane-world theories. Although the dark radiation term 2m/a 4 is 
positive (for m > 0) in the Friedmann equation (16), it appears with negative sign in the 
Raychaudhuri equation (15). Therefore dark radiation increases both the cosmic expansion 
rate and the cosmic deceleration as compared to standard cosmology. A slower expansion 
rate and additional acceleration is gained for negative values of m, but such a possibility is 
not allowed by our present model. 

Would the bulk radiation pushing away the brane and therefore causing cosmic expan- 
sion be able to produce inflation on the brane for any set of numerical data? For the range 
of numerical parameters employed in the preceding Section such an effect does not occur. 
Without radiation, the scale factor shows the typical evolution of a closed universe; cf Fig 2. 
Inflation would mean the occurrence of an exponential expansion, but this is not happening 
by switching on the radiation. Indeed, a glance on the scales of Figs 2 and 8, respectively 
(in the latter we have represented the differences in the scale factor with and without radi- 
ation), shows that radiation is not changing significantly the evolution of the scale factor. 
The deviations of physical quantities evaluated in the radiating case as compared to the 
nonradiating case are typically of the order of 10~ 4 . We have also checked that the radia- 
tion is not changing significantly the evolution of the Hubble-parameter depicted in Fig 3, 
therefore inflation due to this mechanism is excluded. 

Stated more generically and independently of the specific numerical choices in our analy- 
sis, we can pose the question, whether the absorbed radiation could behave like a scalar field 
in a slow-roll regime. Eq. (27) shows that the the absorbed radiation is better interpreted 
as dark dust matter, rather than dark radiation. Indeed, it is a pressureless fluid. Therefore, 
the attempt to describe it in terms of a scalar field would yield ~ 2V and the slow-roll 
condition could not be satisfied. 

As we are considering closed cosmological models with k — 1, it is obvious that we should 
ask whether the radiation pressure would be able to stop the recollapse of the brane-world 
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universe, producing a bounce? Our numerical investigations suggest a negative answer. 
For the whole range of initial brane energy densities that we have checked, the recollapse 
will inevitably occur. The explanation is again the double role of the radiation: it gives 
rise to both a radiation pressure, which is accelerating the motion of brane in the bulk, 
and an increase of brane self-gravity, which is decelerating it. As consequence, the Hubble 
parameter is zero only once during the whole cosmological evolution, when the expansion 
ceases and the recollapse begins. 

Although radiation pressure can in principle accelerate the cosmological expansion, it will 
always encounter the competing effect of the brane becoming heavier as it absorbes more 
and more radiation. In our model radiation pressure acts as dark energy during the whole 
cosmological evolution; however, dark matter (some form of dust) is produced as well - cf 
Eq. (27). This is counterbalancing the accelerating effect. 

Our numerical investigations show, that the increase in the energy density p rad , which 
is due to the absorbed radiation is not significant as compared to the energy density of 
conventional matter. This can be seen from our Fig 9, representing the variation in p due 
to the bulk black hole radiation. The variation is of order 10~ 4 as compared to the energy 
density of the brane in the non-radiating case, shown in Fig 3. It is important to mention 
this point, as in our model we have studied a radiation dominated universe. By contrast, all 
dark matter born from the Hawking radiation absorbed on the brane is dust, cf. Eq. (27). As 
the absorbed energy density merely represents a small fraction of the original brane energy 
density, our assumption of a radiation-dominated universe holds true. 

As we have seen, even the tiny variations of order 1CT 4 in the physical quantities can be 
achieved only at the price of having a huge brane tension A = 10 60 TeV 4 . With the Hawking 
radiation energy density so small, we do not expect drastical changes to occur in the mass 
of the bulk black hole, although theoretically a complete evaporation may occur at tc (see 
the remarks following Eq. (31)). Indeed, the coefficient of (t — t ) in Eq. (31) is very small, 
therefore the mass of the brane will not increase significantly due to the absorbtion of this 
radiation, even seen during the whole cosmological evolution. Accordingly, the bulk black 
hole will not completely evaporate, rather its mass function will decrease only slightly during 
the lifetime of the closed universe. How small the rate of evaporation of the bulk black hole 
is can be better seen from Eq. (38), containing dimensionless variables. There, with our 
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choice M T /M P = 1/10, the coefficient of (t — t ) can be evaluated to give 

*P = ™M x 10 - ~ 7.475 x ID"' . (46) 
V6 3V67T 7 

This means that starting with the initial (dimensionless) mass function of the bulk black 
hole of m L (0) = 0.01, towards the end of the recollapse the mass function decreases to the 
values presented in Table I. We note that he value of the initial brane energy density has 
a certain influence on the final value of the dimensionless mass function, the evaporation 
being maximal for the critical brane energy density. 

TABLE I: Due to Hawking radiation, the dimensionless mass function of the bulk black hole 
slightly decreases during cosmological evolution. Towards the end of the recollapse (at t = 1.115) 
the dimensionless bulk black hole mass function is minimal for the critical value of the initial brane 
energy density. 





p = 100 


po = 520 


po = 2000 


t = 


0.01 


0.01 


0.01 


t = 1.115 


0.0099582386 


0.0099582385 


0.0099582387 



VII. CONCLUDING REMARKS 

In this paper we have considered a scenario which is able to deal with the possibility of 
a Hawking radiation escaping a bulk black hole and we have investigated the consequences 
of such a radiation on the cosmological evolution of a closed k — 1 universe. 

For this purpose first we have developed a formalism suitable for describing the asym- 
metric brane-world containing a single radiating bulk black hole on the left bulk region. In 
the most generic scenario the radiation can be partially absorbed, partially transmitted and 
partially reflected on the brane. We have considered here the case of negligible reflection, 
because we wanted to dispose of an exact solution of the Einstein equations in the bulk. By 
suppressing the reflection, and considering the geometrical optics limit of the radiation, the 
bulk can be described by the VAdS5 metric. 
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The part of radiation absorbed on the brane generates quite interesting effects, which were 
not studied before. These can be related to both dark energy and dark mattter. In order 
to concentrate on these effects, the absorption was chosen as maximal. By suppressing 
transmission, the equations could be partially integrated, cf. Eq. (31). Then we have 
introduced dimensionless quantities in order to perform a complex numerical analysis. 

We have carefully compared the cosmological evolutions of the models both without and 
with radiation and found that the Hawking radiation of the bulk black hole represents merely 
a perturbation of the non-radiating model. This holds true during the entire cosmological 
evolution. 

We have studied these perturbations in detail and shown that in the presence of the 
radiation, there are two competing effects. The absorption on the brane contributes towards 
the self-gravity of the brane, which then tends to recollapse faster. This phenomenon appears 
as continuously emerging dark matter on the brane. By contrast, radiation pressure tends 
to accelerate away the brane from the bulk black hole, producing an equivalent of dark 
energy. For the numerical data employed, we have found the critical value of the initial 
brane energy density for which these two competing effects roughly annihilate each other. 
We have also analysed cosmological evolution for both lighter and heavier branes and found 
that (at early times) the radiation pressure is dominant for light branes, while for heavy 
branes the self-gravity dominates during the whole cosmological evolution. 

As a main result, we have proved that the asymmetry introduced in the model is not 
able to change the recollapse of the universe in the k = 1 case, regardless of whether the 
radiation is switched off or on. Our numerical analysis has shown that the bulk Hawking 
radiation cannot change the final fate of the recollapsing universe either. 
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